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Abstract 


This  report  describes  experimental  and  analytical  investigations  in 
the  area  of  wear  of  materials  due  to  repeated  impact  loading.  A recipro- 
cating impact  wear  test  apparatus  is  utilized  to  provide  controlled, 
repetitive  impacts  between  a material  specimen  and  a counterface.  The 
impulsive  loading  may  be  purely  normal  or  a condition  of  relative  sliding 
may  be  introduced  during  the  impact  event.  A wide  variety  of  materials 
have  been  investigated  in  both  normal  and  sliding  impact  modes.  These 
include  polymers,  composites,  and  several  metals.  The  wear  behavior  of 
these  materials  is  investigated  as  a function  of  peak  impulsive  load, 
relative  sliding  velocity,  and  number  of  impact  load  cycles.  Weight 
loss  and/or  surface  profile  data  are  obtained  as  a function  of  these 
parameters.  Scanning  electron  microscope  studies  of  the  wear  surface 
and  of  sub-surface  sections  suggest  that  delamination  is  the  operative 
mechanism  for  some  materials  but  not  for  others. 

Analytical  investigations  are  directed  toward  determination  of  the 
states  of  stress  in  material  specimens  undergoing  controlled,  repetitive 
impact.  Results  are  derived  from  the  theory  of  elasticity  by  application 
of  the  Papkovich-Neuber  solution  combined  with  the  double,  finite  Fourier 
transformation.  Correlations  between  sub-surface  stress  levels  and  wear 
debris  formation  are  sought. 
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Experimental 

As  noted  earlier  [1],  we  have  studied  a broad  range  of  materials  of 
technological  significance.  These  include  polymers,  composites,  steels, 
and  titanium  and  aluminum  alloys.  The  following  sections  will  summarize 
salient  experimental  results  for  each  category  of  material  noted  above. 

A later  section  will  summarize  results  which  overlap  the  material  domains. 

A.  Polymers 

The  polymeric  materials  tested  were  Teflon  and  Delrin  (polytetra- 
fluoroethylene  and  acetal  resin).  Results  from  these  test  series  have 
been  reported  [2],  It  was  found  that  wear  laminae  were  particularly 
apparent  in  compound  impact,  v here  relative  transverse  sliding  accompanies 
the  normal  impulse.  Increased  relative  sliding  speed  was  found  to 
accelerate  the  process  of  surface  and  subsurface  crack  initiation,  with 
increasing  numbers  of  impact  cycles  tending  to  refine  microcrack  patterns. 

The  relatively  amorphous  nature  of  teflon  as  contrasted  to  the  . 
relatively  crystalline  structure  of  delrin  gave  rise  to  some  concern 
over  the  mechanism  of  crack  initiation  in  teflon.  Based  upon  observa- 
tions of  microvoids  in  commercially  produced  polymeric  materials  [3,4], 
a mechanism  for  crack  initiation  in  relatively  amorphous  polymers  was 
postulated. 

B.  Composites 

The  graphite  epoxy  composites  tested  were  machined  from  a press 
cured  composite  panel  supplied  from  3M  Company  and  from  pultruded  rod 
furnished  by  Hercules,  Inc.  Preliminary  results  on  the  3M  product  have 
been  reported  [5],  and  further  work  with  both  composites  supports  the 
earlier  findings. 
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With  fiber  orientation  parallel  to  the  normal  impact  direction, 
both  composites  are  found  to  possess  excellent  wear  resistance.  In 
particular,  for  compound  impact  loading  where  wear' rates  are  normally 
appreciable  for  metallic  and  polymeric  materials,  the  composites  show 
virtually  no  wear.  For  conditions  of  relatively  high  transverse  sliding 
speeds,  a thin  film  of  epoxy  forms  due  to  interfacial  forces  and  high 
surface  temperatures  arising  in  relative  sliding.  As  seen  in  Figure  1, 
this  film  tends  to  cover  the  ends  of  the  graphite  fibers,  and  acts  as  a 
solid  film  lubricant.  For  the  minimal  wear  that  does  occur,  the  epoxy 
film  appears  to  peel  off  in  thin  wear  sheets  suggestive  of  delamination. 
Figure  2 shows  a composite  view  of  surface  and  subsurface  features,  and 
clearly  demonstrates  the  existence  of  the  epoxy  layer.  Figure  3,  taken 
at  9,000  X,  indicates  the  process  by  which  the  epoxy  resin  flows  from 
within  the  matrix  to  form  the  surface  film. 

C.  Steels 

AF-1410  steel  was  tested  in  two  microstructural  conditions.  This 
steel  (a  modified  HY  130),  is  a weldable,  fracture-resistant,  high 
strength  alloy  developed  for  airframe  structural  applications.  For 
purposes  of  wear  testing,  the  material  was  double  austenitized,  quenched 
at  each  interim,  and  appropriately  aged  to  obtain  the  desired  micro- 
structure. In  one  case,  the  aging  was  conducted  to  obtain  the  maximum 
strength  condition  ("strong"),  while  in  the  other,  the  material  was  aged 
to  yield  the  maximum  toughness  condition  ("tough").  Specimens  from  both 
microstructural  conditions  were  tested  in  compound  impact  at  both  low 
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velocity  (1.4  fps)  and  high  velocity  (12.2  fps).  Peak  normal  impact 
stress  levels  for  these  two  series  were  980  psi  and  6100  psi , respectively. 

Weight  loss  curves  as  a function  of  number  of  cycles  of  impact  for 
the  high  speed  series  are  shown  in  Figure  4.  These  results  demonstrate 
that  "two  body"  wear  is  occurring  for  this  material  in  both  microstructural 
conditions.  The  impact  caps  (counterface  material),  in  fact,  suffer 
considerably  more  wear  than  do  the  specimens.  Further,  one  notes  that 
the  "tough"  microstructure  is  more  wear  resistant  than  the  "strong",  and 
also  produces  less  wear  on  the  counterface. 

Surface  and  subsurface  microscopy  reveals  several  features  of 
importance  for  these  high  velocity  test  series.  Figure  5 shows  a 
representative  "tough"  specimen  surface  following  50,000  cycles  of 
impact  loading.  Higher  magnifications  of  surface  features  give  insight 
into  the  operative  wear  mechanisms.  Figure  6,  for  example,  shows  both 
evidence  of  plowing  (the  relatively  smooth  wear  groove)  and  adjacent 
areas  suggestive  of  delamination  with  characteristic  layering.  Figure 
7,  on  the  other  hand,  is  suggestive  of  delamination  alone,  with  a laminar 
wear  particle  lifting  from  the  surface,  and  about  to  "peel  off". 

Figure  8 conclusively  demonstrates  the  layering  phenomenon,  of  which 
more  will  be  said  later. 

But  surface  features  alone  are  insufficient  to  determine  whether 
the  delamination  mechanism,  as  postulated  in  the  theory  [6],  is  operative. 
Figure  9 shows  a subsurface  section  of  a "strong"  specimen  following 
150,000  loading  cycles.  The  lath  martensite  orientation  indicates  the 
high  degree  of  subsurface  plastic  deformation  which  has  occurred,  yet 
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there  is  no  evidence  of  voids  or  cracks  in  this  region.  The  surface 
layer  itself  is  of  considerable  interest.  One  notes  an  almost  structure- 
less zone  which  is  thought  to  be  a highly  refined  recrystallized  layer. 

There  appears  to  be  some  suggestion  of  cracks  within  this  layer,  and  on 
the  interface  between  it  and  the  severely  plastically  deformed  substrate. 

More  will  be  said  later  concerning  these  subsurface  zones. 

The  low  velocity  test  series  for  the  AF -1410  steels  yielded  several 
results  of  significance.  First,  "two  body"  wear  was  observed,  in  which 
both  specimen  ( AF- 1410)  and  counterface  (impact  cap)  experienced  considerable 
wear.  The  counterface  material  is  a 17-4  PH  stainless  steel,  and  is 
described  in  more  detail  below.  From  Figure  10,  one  notes  that  the  1410 
steels  produced  several  times  more  counterface  wear  than  did  the  titanium 
alloys,  and  orders  of  magnitude  more  wear  than  did  the  aluminum  alloys. 
Indeed,  for  the  1410  steels,  the  counterface  wear  well  exceeded  the 
specimen  wear,  suggesting  that  the  1410  steels  be  considered  as  good 
wear  resistant  materials. 

In  addition,  as  can  be  seen  in  Figure  11,  the  1410  specimens  demon- 
strate significantly  different  wear  rates.  The  1410  "strong"  shows  an 
initially  high  wear  rate  up  to  50,000  load  cycles,  with  a very  modest, 
steady  state  wear  rate  up  to  160,000  cycles.  In  contrast,  the  1410 
"tough"  possesses  very  modest  wear  up  to  75,000  cycles,  with  an  increasing 
wear  rate  up  to  160,000  impact  load  cycles.  At  this  outer  limit  for 
this  test  series,  the  actual  specimen  wear  overlaps  for  these  two  micro- 
structures at  160,000  cycles.  From  both  Figures  10  and  11  it  is  clear 
that  microstructural  variables  play  an  important  role  in  impact  wear. 

The  counterface  material  utilized  in  all  tests  with  all  specimen 


materials  is  17-4  PH.  This  material  can  be  broadly  classified  as  a 
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precipitation  hardening  martensitic  stainless  steel.  Following  the 
finding  that  this  material  was  subject  to  considerable  wear  with  particular 
specimen  materials  (most  notable  the  1410  steels,  but  also  the  titanium 
alloys),  surface  and  subsurface  microscopy  were  conducted  on  selected 
impact  caps.  Figure  12  shows  a surface  micrograph  of  such  a counterface 
when  run  against  1410  "tough"  in  high  speed  sliding  impact  for  150,000 
load  cycles.  In  this  photo,  the  base  17-4  PH  material  can  be  seen  as 
the  bottom  surface  of  a crater  from  which  a wear  particle  has  been 
plucked.  The  surrounding  zones  indicate  plowing  striations  and  a large 
laminar  particle  layered  over  this  plowed  region.  These  features,  taken 
together,  suggest  that  some  material  transfer  is  occurring  between 
specimen  (1410)  and  counterface  (17-4  PH). 

Subsurface  microscopy  on  impact  caps  is  in  progress.  Preliminary 
findings  from  a low  velocity  1410  test  show  evidence  of  a thin  surface 
layer  on  the  impact  caps,  but  there  is  no  evidence  of  void  or  crack 
formation,  nor  even  any  significant  subsurface  plastic  deformation  in 
the  impact  cap  substrate. 

D.  Titanium  Alloys 

Two  titanium  alloys  were  tested  in  both  high  and  low  velocity 
compound  impact  series.  The  alloys  are  IMI  685  and  RMI  5522S;  the 
compositions  of  these  alloys  have  been  reported  earlier  [1],  as  have 
some  preliminary  results  [7].  Both  of  these  alloys  have  been  developed 
for  high  temperature  service,  in  which  strength  and  creep  resistance  are 
important. 

For  purposes  of  wear  testing,  the  materials  were  processed  as 
follows.  The  RMI  5522S  was  received  as  1/2  inch  diameter  swaged  rod, 
and  was  not  heat  treated  in  order  that  a microstructure  consisting  of 
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large  amounts  (about  90%)  of  continuous  elongated  primary  alpha  might  be 
tested  for  impact  wear  resistance.  While  this  microstructure  is  far 
from  optimum  in  terms  of  high  temperature  creep  resistance,  it  represented 
an  unknown  in  terms  of  wear  resistance. 

The  I MI  685  alloy,  on  the  other  hand,  was  processed  to  its  re- 
commended microstructural  condition,  consisting  of  nearly  100%  trans- 
formed beta.  This  material  was  received  as  1 inch  diameter  rod,  and  was 
subsequently  swaged  to  1/2  inch  diameter;  it  was  then  beta  annealed, 
quenched,  and  heat  treated  to  produce  the  typical  "basketweave"  micro- 
structure. 

In  high  velocity  compound  impact  testing,  the  IMI  685  alloy  showed 
only  a minor  velocity  dependence,  while  the  RMI  5522S  demonstrated  a 
significant  dependence  as  seen  in  Figure  13.  The  peak  normal  stress  for 
this  test  was  1500  psi,  and  it  is  found  that  more  wear  occurs  at  the 
lower  transverse  sliding  velocity  (10.5  fps)  than  at  the  higher  (17.5 
fps).  While  this  finding  conflicts  with  results  obtained  in  ballistic 
impact  wear  studies  [8],  it  agrees  with  recent  results  in  pure  sliding 
wear  (no  impact)  obtained  with  commercially  pure  titanium  at  high  sliding 
speeds  [9].  In  this  latter  study,  wear  rates  were  found  to  exhibit  both 
a minimum  and  a maximum  within  a range  of  sliding  velocities.  Between 
1.3  and  3.3  fps,  the  wear  rate  decreased;  between  3.3  and  16.4  fps,  the 
wear  rate  increased;  and  as  the  sliding  speed  was  increased  from  16.4  to 
32.8  fps,  the  wear  rate  again  decreased. 

It  must  be  appreciated  that  in  this  regime  of  high  speed  sliding, 
significantly  high  surface  temperatures  are  generated,  both  in  sliding 
and  compound  impact  wear.  Indeed,  both  in  [9]  and  in  our  own  work  with 
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titanium  alloys,  sparks  are  observed  during  testing.  These  high  surface 
temperatures  are  thought  to  lead  to  a local  softening  of  the  specimen 
surface,  with  a concomitant  reduction  in  frictional  force  at  the  interface. 
Thus,  at  the  higher  sliding  velocity,  reduced  friction  is  believed  to 
account  for  the  lesser  degree  of  wear  observed.  This  reasoning,  however, 
cannot  be  extrapolated  to  a low  velocity  test  series,  where  thermal 
effects  are  less  significant. 

Figure  14  is  representative  of  the  surface  appearance  of  both 
titanium  alloys  following  high  speed  sliding  impact.  While  the  flat 
sheets  are  suggestive  of  delamination,  there  is  also  the  plowing  groove 
network  superposed.  It  is  possible,  indeed,  that  the  flat  particles  are 
those  which  have  been  back-transferred  from  the  counterface,  and  have 
been  pounded  to  shape  during  impact. 

More  interesting  are  the  subsurface  photomicrographs  from  the  high 
speed  sliding  tests.  Figure  15  is  at  relatively  low  magnification  and 
shows  the  specimen  sectioned  parallel  to  the  sliding  direction.  One 
notes  a leading  and  a trailing  edge  of  material  which  indicate  the  high 
degree  of  plastic  deformation  which  occurs.  Figure  16,  a higher  magnifi- 
cation view  of  the  same  material  (IMI  685),  shows  a plastically  deformed 
surface  layer  which  contains  a swatch  of  material  which  is  almost  broken 
free.  Figure  17  is  a subsurface  section  from  RMI  5522S  and  demonstrates 
the  three  zones  typically  found  in  such  examinations.  The  deep  substrate 
consists  of  the  unaffected  microstructure  (zone  1).  The  intermediate 
layer  is  plastically  deformed,  and  this  deformation  is  typically  observable 
in  the  microstructure  (zone  2).  Often  a refinement  in  this  microstructure 
can  be  observed  as  the  near-surface  layer  (zone  3)  is  approached. 

Zone  3 is  thought  to  be  a severely  plastically  deformed,  recrystallized 
layer.  In  this  case,  obvious  voids  and  cracks  are  visible  in  zone  3. 
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In  the  low  velocity  test  series  on  the  titanium  alloys,  transverse 
sliding  speed  was  1.4  fps  and  peak  normal  stress  level  was  980  psi.  As  j 

noted  earlier,  both  titanium  alloys  gave  rise  to  counterface  wear,  but  I 

not  to  the  degree  that  the  AF-1410  steels  did  (Figure  10).  Also  as 
noted  earlier  (Figure  11),  the  titanium  alloys  exhibited  considerable 
wear  in  the  low  velocity  test  series.  Indeed,  these  materials  (and 
aluminum  alloy  2011-T3)  were  subject  to  significantly  higher  wear  rates 
than  2124  aluminum,  the  AF-1410  steels,  and  the  graphite  epoxy  composites 
(not  shown  in  Figure  11  as  the  wear  curve  is  indistinguishable  from  the 
abcissa). 

Surface  microscopy  for  the  low  speed  sliding  impact  test  series 
reveals  little  difference  from  results  observed  in  high  speed  compound 
impact.  Figure  18  is  representative  for  low  speed,  and  should  be  compared 
to  Figure  14  for  high  speed.  Subsurface  results  are  also  similar  between 
low  and  high  speed  sliding,  as  seen  in  Figures  19  and  20.  Figure  19 
shows  the  RMI  5522S  alloy,  where  again  the  three  characteristic  zones 
can  be  observed,  while  Figure  20  shows  the  IMI  685  alloy  at  high  magnification. 

This  material  showed  virtually  no  observable  subsurface  features  from 
the  low  speed  test  series. 

Taken  together,  and  noting  significant  transfer  films  and  counter- 
face wear,  it  appears  that  not  only  delamination  but  also  adhesive  and 
plowing  mechanisms  are  operative  with  these  titanium  alloys  over  the 
range  of  compound  impact  parameters  explored. 

E.  Aluminum  Alloys 

Two  aluminum  alloys  have  been  tested  in  both  low  and  high  speed 
compound  impact.  The  alloys  are  2124  and  2011-T3.  The  2124  is  a wrought 
product  consisting  of  aluminum-copper  solid  solution  with  a dispersion 
of  intermetal  lies  contained  within  the  relatively  soft,  tough  matrix. 
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This  alloy  also  possesses  very  good  fracture  toughness,  particularly  in 
the  short  transverse  direction  from  which  impact  specimens  were  machined. 

An  alternate  designation  for  2124  is  2024-T851. 

The  2011-T3  alloy  is  a extruded  product  which  consists  of  aluminum- 
copper  solid  solution  matrix  in  which  there  are  dispersed  particles  of 
lead  and  bismuth.  This  alloy  is  "free  machining",  as  the  soft,  low 
melting  point  lead  and  bismuth  constituents  bond  poorly  to  the  matrix 
and  act  as  incipient  voids  or  crack  initiators. 

High  speed  sliding  tests  with  these  alloys  were  conducted  at  a 
variety  of  velocities  and  peak  stress  levels.  Sliding  velocities  ranged 
from  8.7  to  19.2  fps;  peak  stress  levels  ranged  from  610  to  3260  psi. 

Wear  rates  were  quite  high  for  the  2011 -T3  and  moderate  for  the  2124. 

Surface  SEM  photomicrographs  for  the  two  alloys  are  virtually 
interchangeable.  Accordingly,  only  those  from  the  2124  alloy  are  presented. 
Figure  21  shows  a relatively  low  magnification  view  of  the  surface  of  a 
specimen  run  at  19.2  fps  for  50,000  load  cycles.  On  this  level,  the 
wear  mechanism  appears  to  be  abrasive  as  evidenced  by  the  deep  plowing 
grooves.  At  increased  magnification,  however,  as  seen  in  Figure  22, 
there  is  evidence  of  layering  and  of  flat,  lamellar  particles.  Figures  23 
and  24  show  such  features  at  higher  magnifications. 

Subsurface  micrographs  for  this  2124  alloy  are  shown  in  Figures  25 
and  26.  Figure  25  is  typical  for  this  material,  and  subsurface  features 
are  difficult  to  distinguish  in  the  scanning  electron  microscope.  There 
is  little  evidence  of  a region  of  large  plastic  deformation  (zone  2), 
although  some  flow  patterns  can  be  distinguished.  The  near-surface 
region  (zone  3)  is  again  clearly  demarked,  and  Figure  25  suggests  that 
it  is  separated  from  zone  2 by  a propagating  crack.  Figure  26  shows  a 
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dark  field  image  from  the  light  microscope.  This  photo  clearly  indicates 
the  gross  plastic  deformation  which  has  occurred  in  producing  the  curled 
"wear  lip"  on  the  trailing  edge  of  the  specimen. 

But  while  there  is  an  apparent  absence  of  subsurface  features  in 
the  2124  alloy,  the  2011-T3  contains  an  abundance.  Figure  27  is  a photo 
taken  at  400X  from  a specimen  run  at  a sliding  velocity  of  17.5  fps  and 
peak  normal  impact  stress  of  2450  psi  at  250,000  load  cycles.  This 
subsurface  section  is  normal  to  the  sliding  direction,  so  that  sliding 
occurred  into  the  photo.  The  curvature  of  the  surface  is  due  to  the 
more  macroscopic  undulations  observable  as  "plowing  grooves"  at  lower 
magnifications.  One  notes  the  appearance  of  numerous  voids  and  cracks 
in  the  near-surface  zone. 

Figures  28  and  29  show  these  features  in  more  detail.  Figure  28  is 
a section  taken  parallel  to  the  sliding  direction,  where  sliding  occurred 
from  right  to  left.  One  notes  numerous  voids  and  cracks,  suggestive  of 
subsurface  features  found  in  support  of  the  delamination  theory  of  wear. 
Figure  29  is  again  normal  to  the  sliding  direction,  and  here  it  appears 
that  one  near-surface  zone  has  been  deposited  over  another. 

Figure  30  shows  another  section  parallel  to  the  sliding  direction, 
and  again  cracks  run  essentially  parallel  to  the  surface.  Figure  31  is 
included  to  demonstrate  the  non-uniformity  of  subsurface  features.  In 
this  micrograph,  there  is  a comparative  absence  of  the  marked  void  and 
crack  features  noted  above.  However,  one  notes  a distinct  near-surface 
layer  which  appears  to  be  separating  from  the  plastically  deformed 
zone  2.  Figure  32  shows  a parallel  section  from  a specimen  which  was 
run  for  only  10,000  cycles.  Thus,  one  finds  significant  subsurface  void 
and  crack  formation  in  the  2011-T3  alloy  after  a relatively  few  impact 
load  cycles. 
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The  low  speed  sliding  tests  with  the  aluminum  alloys  result  in 
similar  findings.  As  shown  in  Figure  10,  virtually  no  wear  of  the 
counterface  (impact  caps)  occurs  with  either  alloy.  And,  as  depicted  in 
Figure  11,  the  2124  alloy  is  considerably  more  impact  wear  resistant 
than  is  the  2011-T3. 

All  of  these  results  are,  those  which  one  would  expect  from  a consideration 
of  the  microstructures  of  the  two  alloys.  While  the  2124  is  formulated 
for  high  fracture  toughness,  the  2011-T3  is  designed  to  facilitate  crack 
formation.  The  subsurface  features  observed  correlate  with  these  differences 
in  microstructure,  and  the  observed  differences  in  wear  rates  likewise 
correlate. 

F.  Reproducibility 

A series  of  tests  were  undertaken  to  determine  the  scatter  in 
weight  loss  data  in  compound  impact.  Specimens  of  5/32"  diameter  201 1 - 
T3  alloy  were  run  at  peak  impulsive  stress  of  1566  psi  with  transverse 
velocity  of  8.7  fps.  Results  were  obtained  from  three  to  six  specimens 
at  each  of  the  following  number  of  impact  load  cycles:  10,000;  50,000; 

100,000;  250,000.  A typical  example  of  the  spread  of  data  is  at  N = 

100,000  load  cycles.  Here,  the  mean  of  six  data  points  is  0.01025  g. , 
and  the  standard  deviation  is  0.00260  g. 

G.  Velocity  Dependence 

Figure  33  presents  weight  loss  data  for  the  2011-T3  alloy  at  two 
transverse  sliding  velocities,  v = 5.5  fps  and  v = 8.7  fps.  As  above, 
specimen  diameter  is  5/32"  and  peak  impulsive  stress  is  1566  psi. 

The  observations  are  that  more  wear  occurs  at  the  lower  relative 


sliding  speed,  and  that  weight  losses  increase  with  increasing  numbers 
of  impact  cycles.  These  results  correspond  to  those  obtained  with 
titanium  alloys. 
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H.  Stress  Dependence 

A special  test  series  was  undertaken  to  determine  the  effect  of 
peak  impact  stress  on  the  wear  rate  of  aluminum  2011-T3.  In  this  series, 
the  peak  impact  force  was  maintained  at  a constant  value  of  75  pounds 
and  the  relative  sliding  velocity  was  likewise  maintained  constant  for 
all  tests  at  17.5  fps.  By  changing  specimen  diameters,  different  peak 
impact  stress  levels  were  imposed.  Figure  34  shows  some  results  from 
this  series,  and  the  stress  dependency  is  obvious:  wear  increases  with 
increasing  stress. 

The  data  from  this  test  series  are  plotted  as  change  in  height  of 
the  specimen  as  a function  of  number  of  load  cycles.  This  is  done  for 
the  following  reason.  If  one  considers  weight  loss  (or  volume  change), 
the  larger  diameter  specimens  would  lose  more  weight  (or  volume)  per 
unit  thickness  worn  away  than  would  the  smaller.  Thus,  the  height 
change  is  more  indicative  of  the  amount  of  wear,  and  is  the  quantity  of 
practical  importance  as  well.  Weight  loss  comparisons  are  significant 
when  one  is  comparing  similar  materials  of  the  same  geometrical  shape 
(such  as  two  titanium  alloys  of  the  same  diameter).  But  when  one  attempts 
to  draw  comparisons  between  materials,  or  between  specimens  of  different 
diameter,  height  change  is  the  significant  variable. 

I.  Miscellaneous  Tests 

In  addition  to  test  results  noted  above,  other  avenues  have  been 
and  are  being  pursued.  These  are  briefly  noted  below. 

Profilometry  has  been  employed  in  two  areas.  One  such  investigation 
concerned  the  determination  of  surface  roughness  of  specimens  in  various 
stages  of  wear.  Another  utilization  of  the  profilometer  is  in  obtaining 
wear  data  for  the  graphite  epoxy  composites.  These  materials  are  exceedingly 
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difficult  on  which  to  perform  weight  less  measurements.  This  is  due  to 
the  fact  that  they  absorb  water  from  the  atmosphere.  Accordingly,  in 
some  instances  we  have  utilized  round  nosed  (spherically  ended)  specimens 
which  develop  a flat  upon  wear  testing.  Profilometer  measurements  allow 
calculation  of  the  worn  volume,  from  which  comparisons  between  the 
composites  and  other  materials  may  be  drawn. 

Other  areas  of  experimentation  include  the  measurement  of  specimen 
surface  temperatures.  This  is  being  accomplished  by  means  of  thermocouples 
embedded  in  a test  material  just  under  the  wear  surface.  As  the  wear 
process  progresses,  the  thermocouple  junction  nears  the  actual  wear 
surface.  Finally,  the  junction  becomes  a part  of  the  surface,  and  the 
test  self-destructs.  In  the  process,  however,  one  obtains  data  on 
temperatures  in  and  on  the  wear  surface.  It  is  felt  that  such  measurements 
are  important,  and  due  to  the  wide  variations  reported  in  the  literature, 
analytical  methods  are  thought  to  be  unreliable. 

Also  in  progress  are  test  series  on  two  additional  materials. 

These  include  precipitation-hardening  copper-chromium  alloys  aged  for 
different  periods  of  time  to  produce  differing  second-phase  particle 
structures.  Similar  tests  will  be  done  with  spherodized  1045  steel. 

These  materials  will  contain  varying  volume  fractions  of  second-phase 
particles,  and  the  particles  themselves  will  vary  in  size  and  spacing. 
Quantitative  metallography  will  be  conducted  in  concert  with  these  test 
series  in  order  to  attempt  to  obtain  quantitative  measures  of  void  and 
crack  formation. 

Summary  of  Experimental  Results 

From  the  results  obtained  to  date,  we  can  note  several  important 
generalizations.  First  among  these  is  the  importance  of  external  variables 
in  the  impact  wear  process,  including  peak  impulsive  stress  levels  and 
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relative  transverse  sliding  velocity.  Secondly,  it  is  clear  that  different 
materials  possess  widely  varying  impact  wear  resistance  which  would  not 
be  inferred  from  simple  comparisons  of  hardness,  creep  resistance,  or 
the  like.  Within  a given  category  of  materials,  microstructural  variables 
are  most  significant.  Thus,  one  should  in  principle  be  able  to  work 
toward  a material  structure  most  favorable  for  impact  wear  resistance, 
and  perhaps  by  extension,  to  wear  resistance  in  other  forms  of  loading. 

With  metallic  materials,  subsurface  phenomena  appear  to  be  significant. 
In  general,  one  observes  three  characteristic  zones  beneath  the  wear 
track.  The  zone  furthest  from  the  wear  surface  ranges  from  unaffected 
microstructure  toward  observably  deformed  "grains".  The  intermediate 
zone  consists  of  material  which  has  been  subject  to  exceedingly  high 
plastic  deformation,  but  the  basic  "grain"  structure  is  still  observable, 
although  deforming  and  refining  into  flow  lines.  The  near-surface  zone 
is  still  more  plastically  deformed,  and  probably  consists  of  material 
both  from  the  specimen  and  that  deposited  from  the  counterface.  It 
appears  to  be  highly  refined  and  homogenized,  and  preliminary  TEM  results 
of  our  own  indicate  it  to  be  crystalline. 

Investigations  by  other  workers  have  recently  been  conducted  on 
such  near-surface  layers  arising  in  simple  sliding  (no  impact)  wear 
[10-12].  Both  Bill  [10]  and  Van  Dijck  [11]  find  this  zone  to  be  crystalline, 
as  distinct  from  the  amorphous  Beil  by  layer  which  results  from  polishing 
and  is  often  cited  in  the  literature.  In  particular.  Van  Dijck  finds 
the  smallest  grains  ever  measured,  ranging  in  thickness  from  200  to  600 

O 

A.  Moreover,  the  plastic  deformations  in  this  layer  are  likewise  the 
largest  ever  measured,  ranging  up  to  strains  of  230. 

Clearly,  the  formation  and  properties  of  these  zones  occurring  in 


J 


impact  wear  need  further  study. 
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Analytical 

During  the  second  phase  of  our  work,  a more  realistic  analytical 
description  of  the  wear  test  and  a better  physical  interpretation  of  the 
analysis  was  achieved.  As  noted  earlier  [1],  we  are  dealing  with  "controlled 
impact  differing  from  various  forms  of  free  (or  ballistic)  impact  considered 
by  others.  In  essence,  our  impact  wear  testing  machine  forces  the 
specimen  to  strike  the  transducer,  which  may  be  stationary  or  rotating 
transverse  to  the  direction  of  normal  approach.  This  forced  contact  is 
"controlled"  in  the  sense  that  the  specimen  is  driven  forward  by  means 
of  the  spatial  mechanism  of  the  impact  wear  machine  [13].  The  dynamic 
characteristics  of  the  spatial  mechanism,  which  vary  depending  upon 
several  adjustments  made  for  specific  materials  and  specific  test  conditions, 
serve  to  interact  with  the  dynamic  response  characteristics  of  the 
transducer  counterface  and  its  mountings.  In  short,  the  impact  is 
"controlled"  as  compared  to  the  "free"  impact  that  would  be  characterized 
by  a ballistic  device,  or  a free-fall  device. 

Accordingly,  the  analysis  is  complicated  by  the  fact  that  the 
motion  of  the  far  end  of  the  specimen  is  not  known  exactly,  due  to  the 
dynamic  response  characteristics  of  the  overall  system,  and  the  existence 
of  tolerances  between  moving  components.  Hence,  the  unknown  motion 
characteristics  of  the  far  end  of  the  specimen  must  be  determined.  This 
is  possible,  by  means  of  a rather  complicated  procedure  developed  by 
Solecki  during  the  first  phase  of  our  work.  Thus,  the  initial  analysis 
dealt  with  the  impact  of  a homogeneous  one-dimensional,  finite  elastic 
rod  [14]  upon  an  elastic  half-space  (Figure  35).  In  this  work,  it  was 
assumed  that  the  far  end  of  the  rod  was  subjected  to  a prescribed  motion 
due  to  a controlled  impact.  The  friction  between  the  rod  and  the  half- 
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space  was  disregarded  in  order  to  avoid  dynamic  bending  of  the  rod  in 
addition  to  the  propagation  of  longitudinal  stress  waves.  The  objective 
of  this  preliminary  work  was  to  investigate  procedures  useful  when 
considering  the  impact  of  a 3-dimensional  body  and  to  acquire  a better 
understanding  of  the  mechanical  behavior  of  a system  under  conditions  of 
controlled  impact.  This  work  is  now  in  a final  stage  leading  to  numerical 
results,  some  of  which  have  been  reported  earlier  [1].  The  main  points 
of  the  theory,  an  explanation  of  adopted  techniques,  and  the  derivation 
of  analytical  relations  are  briefly  described  in  Appendix  2. 

It  was  found  [see  Appendix  2]  that  the  displacement  of  an  arbitrary 
point  of  the  rod,  u(x,t),  is  given  by  the  expression 

u(x,t)  = V A(x,t)/p7E  + (f[t-x/p7E]  - f[t-(2L-x)^p7E]+...} 

+ {w[t-(L-x)/p/E]  - w[t- (L+w)/p/E]+. . . } 

where  f is  the  "controlling  function";  w is  the  displacement  of  the 
striking  end;  /E/p"  is  the  speed  of  elastic  wave  propagation;  V is  the 
velocity  of  the  rod  at  the  instant  of  the  impact  and  A(x,t)  is  a bi-  or 
trilinear  function  of  x and  t.  An  inspection  of  this  expression  reveals 
qualitatively  the  nature  of  the  state  of  stress  during  the  controlled 
impact.  The  first  term  on  the  right  hand  side,  proportional  to  the 
velocity  of  the  impact,  is  responsible  for  generating  two  rectangular 
stress  waves:  a compression  wave  propagating  from  the  striking  end  and 
a tension  wave  propagating  from  the  far  end.  These  waves  are  reflected 
from  the  opposite  ends  (changing  sign  at  the  far  end)  and  continue 
propagating  until  the  separation  time.  This  results  in  periodically 
oscillating  stress  at  any  cross-section  of  the  rod.  The  second  element 
in  this  expression,  represented  by  the  collection  of  terms  in  the  first 
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brackets,  describes  a finite  series  of  "signals”  generated  at  the  far 
end  by  the  controlling  function,  f(t).  The  effect  of  this  is  a train  of 
stress  waves  transmitted  continuously  by  the  far  end  until  the  time  of 
separation  occurs.  If  the  duration,  t^ , of  the  contact  is  large  compared 
to  the  traverse  time,  T,  of  the  elastic  wave,  then  the  maximum  stress 
resulting  from  this  second  element  differs  only  slightly  from  the  static 
stress.  The  last  element  (collection  of  terms  in  the  second  brackets) 
is  similar  in  its  form  to  the  second  term  and  describes  the  effect  of 
the  elastic  flexibility  of  the  support.  Thus,  a sequence  of  stress 
waves  originating  from  the  striking  end  tends  to  decrease  the  pressure 
(or  tension)  generated  by  the  controlled  motion  of  the  far  end.  It  is 
worth  emphasizing  that  in  the  case  of  the  controlled  impact,  the  velocity, 

V,  is  not  an  independent  quantity  but  is  rather  the  time  derivative  of 
the  function  f(t)  at  the  instant  of  the  impact.  Hence  the  first  two 
elements  are  closely  related  and  a simplified  analysis  indicates  that 
the  importance  of  the  first  "impact"  term  is  amplified  relative  to  the 
second  term  as  the  duration  time  of  the  contact  decreases. 

The  three-dimensional,  cylindrical  model  of  the  striking  specimen 
will  have,  of  course,  a different  stress  distribution  than  that  in  the 
rod.  This  is  due,  in  part,  to  the  presence  of  shear  waves,  reflection 
of  waves  from  the  curved  surfaces  of  the  cylinder,  dispersion,  etc. 

Nonetheless,  the  basic  character  of  the  physical  process  will  be  quali- 
tatively the  same  as  that  described  above. 

It  is  also  expected  that  a quasi-static  counterpart  of  this  three- 
dimensional  model  will  have  a similar  stress  distribution  at  the  striking 
end.  Thus,  during  the  second  phase  of  our  work,  our  primary  effort  was 
concentrated  in  the  following  quasi-static  problem.  A finite,  homogeneous, 
isotropic,  elastic  flat-ended  circular  cylinder  is  subjected  to  a given 


- 18  - 


displacement  at  the  upper  end  (constant  vertical  component  of  the  displacement; 
the  remaining  components  equal  to  zero).  The  lower  end  impacts  upon  a 
moving,  rough,  rigid  half-space.  The  roughness  of  the  half-space  generates 
tangential  stresses  due  to  friction,  while  the  vertical  component  of  the 
displacement  eouals  zero  (Figure  36).  Since  this  problem  is  not  axially 
symmetric,  it's  solution  involves  considerable  difficulties.  The  Papkovich- 
Neuber  functions  were  used  (for  a detailed  explanation  of  the  technique, 
see  Appendix  1)  to  achieve  partial  uncoupling  of  [Javier's  equations  of 
equilibrium  in  a polar  coordinate  system.  This  was  followed  by  double 
finite  Fourier  transformation  with  respect  to  e and  z resulting  in 
complete  uncoupling  and  in  solution  of  the  transformed  differential 
equations.  Application  of  the  boundary  conditions  at  the  stress-free 
curved  surface  of  the  cylinder  and  of  the  conditions  at  its  flat  ends 
leads,  after  numerous,  extremely  cumbersome  rearrangements  to  a system 
of  Fredholm's  integral  equations  of  the  first  kind. 

The  resulting  stress  distributions  are  expected  to  correlate  with 
experimental  observations  of  sub-surface  plastic  deformation  and  void 
and  crack  formation.  For  at  least  some  materials,  these  phenomena  are 
thought  to  be  intimately  related  to  the  formation  of  wear  debris. 
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Appendix  1.  The  effect  of  friction  on  the  stress  distribution  in  a finite 
elastic  cylinder. 

1.  Introduction 

In  the  present  paper  the  effect  of  friction  on  the  stress  distribution  in 
a finite,  isotropic,  homogeneous  elastic  circular  cylinder  is  analysed.  It  is 
assumed  that  the  curved  surface  of  the  cylinder  is  stress-free.  The  components 
of  the  displacement  vector  are  preassigned  at  one  flat  end  of  the  cylinder 
while  at  the  other  end  the  axial  component  of  the  displacement  vector  is  zero. 
Also  at  this  end  the  shear  component  of  the  stress  tensor  is  either  zero  (in 
one  direction)  or  proportional  to  the  normal  component  of  the  stress  tensor 
(in  the  perpendicular  direction).  The  problem  is  thus  not  axially  symmetric. 
The  Papkovich-Neuber  functions  are  used  in  order  to  uncouple  the  field  equa- 
tions. The  dependence  of  the  unknown  functions  on  the  angular  and  axial  vari- 
ables is  then  eliminated  by  means  of  finite  Fourier  integral  transformations. 
Application  of  boundary  conditions  at  the  flat  ends  leads  to  a coupled  system 
of  Fredholm  integral  equations  of  the  first  kind  where  the  unknowns  are  the 
shear  stresses  at  those  ends.  Numerical  techniques  will  be  employed  to  solve 
this  system. 

2.  Solution  of^  the  field  equations 

The  differential  equations  governing  the  displacements  of  a homogenous, 
isotropic  elastic,  cylinder  take  the  following  form  in  the  cylindrical  coor- 
dinate system  (see  [1]  p.  236). 
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where  uf,  u0,  are  the  physical  components  of  the  displacement  vector  V2 
is  the  Laplacian 

,2,  ii  + ii_  + 1 i!  + i!  m 

0r2  ' Sr  r2  062  022 

and  the  dilatation  A is  given  by 
3u  n 3u  3u 

* = ^ + 7 (-39  + V + IF  (3) 


The  partial  uncoupling  of  the  above  system  succeeds  with  the  help  of  Papkovich- 
Neuber  substitution  [2] 


u = l (4  + 1/2  jr-Y)  - 4(.l-v)J£ 


It  is  assumed  without  loss  of  generality  ([2]  p.  331)  that  the  axial  component 
Vz  of  the  vector  potential  ¥ equals  zero  and  therefore  eq  (4)  written  in 
component  form  becomes 
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The  substitution  of  expression  (5)  into  the  Navier  equations  (1)  leads  to  their 
partial  uncoupling 
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The  solutions  of  eq.  6 are  subjected  to  the  following  boundary  conditions  on 
the  curved  surface  r = R: 
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and  the  following  ones  at  z=l: 
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At  the  end  z = 0 the  conditions  o = 0,  a = ko  become 

zy  xz  zz 


o sin  0 + a.  cos  0=0 
rz  0z 


a cos  0 - a.  sin  0 = ka 
rz  0z  zz 


wherefrom,  after  elimination  of  a , one  obtains 
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Therefore  the  conditions  at  z = 0 become  finally 
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The  variable  9 is  eliminated  first  by  applying  finite  Fourier  sine  or  cosine 
transformations : 

2ir 

V = V (r,n,z)  = J'i  (r,0,z)  sin  n0d0 
r r o r 


?0  = ^(r.n.z)  = Jvq  (r,0,z)  cos  n0d0 


<)>=<!>  (r,n,z)  = /<^( r,0,z)  sin  n0d0 


Applying  Fourier  sine  transformation  to  (6.1)  and  (6.3)  and  Fourier  cosine 
transformation  to  (6.2)  results  in  the  following  system  of  equations: 
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where 
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It  can  be  also  shown  that  the  condition  u = const  at  z = const  implies  that 
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at  z = const  also  following  relations  hold 
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The  boundary  conditions  (7.1)  and  (7.2)  are  subjected  to  Fourier  sine  trans- 
formation with  respect  to  9 while  (7.3)  is  subjected  to  Fourier  cosine 
transformation  yielding  at  r = R: 
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Fourier  sine  transforms  with  respect  to  0 of  (8.1)  and  (8.3)  and  Fourier 
cosine  transform  of  (8.2)  have  the  form: 
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Finally  0 is  eliminated  from  the  boundary  conditions  (8.4),  (8.5)  by  means 
of  Fourier  sine  transformation: 
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where 


1 for  n ■ 0 


1/2  for  n | 0 


and  where  the  following  notation  has  been  adopted 
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Elimination  of  the  variable  z follows  by  means  of  finite  Fourier  cosine  transfor- 


mation. 


We  define: 


(r,n,m)  = / ? (r,n,z)  cos  X z dz 
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The  transformed  field  equations  (10)  become  now 
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Also  the  boundary  conditions  are  transformed  again  to  eliminate  z.  One  obtains 
at  r ■ R: 
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Note  that  the  conditions 
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resulting  from  the  absence  of  the  shearing  stress  on  the  curved  surface  have  been 


used  here. 


Eqs.  (17.1)  and  (17.2)  are  uncoupled  by  either  adding  or  substracting  (a 
method  applied  by  Achenbach  [1]  p.  239).  When  the  resulting  equations  are  solved 
one  gets  for  m + 0: 
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while  for  m - 0 the  differential  equations  reduce  to  Euler's  equations  nnd  their 
solutions  have  the  following  forms : 
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where  ¥ and  other  functions  with  a subscript  ''p'1  are  the  particular  solution 
defined  later. 

Using  the  boundary  conditions  eq.  (19)  and  denoting 
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we  obtain  after  numerous  manipulations  the  following  expressions  for  the  con- 
stants C.  , C_  , C.  , H,  , H_  , and  H_  : 
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_ I4v2  + 4v3)y2n2  + 8 (1-v)  (3-12v  + 8x>2)n2  + y®  - 2(l-v)(3-12v 
in  m 

+ 8^>0  Vl  Vl  1 
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+ 2 -5T  x;«  S (ln'’  - * V - 20-'-)  (3-12'-  + I2.! 

+ [n2  + 8v(l-v)n  + y2  - 2(l-v)(3-12v  + 8v2)]y2  I2 

m m n+1 

+ 2[-2n*  - 3y~n2  + 4 (1-v) (3-4v) n2  - y^  + 2(l-v)(3-12v  + 8v2)y2] 
m m m 


1 I 1 . i } 

n-1  n+1 


+ 2^  X - { [ (-l+20v  - 16v2)n2  - 8v(l-v) (3-2v)n  - (l-8v  + 8v2) 

DR2  mn  n 

- y2  + 2(l-v)(3-12v  + 8v2)]y2  I2  , + [(-l+20v  - 16v2)n2  + 8v(l-v)(3 

m m n-1 

-2v)n  + (l-8v  + 8v2)y2  + 2(l-v)(3-12v  + 8v2)b2  I2+1  + 2[2(l-12v 

+ 8v2)n4  + (3-20v  + 16v2)y2n2  - 4(l-v)(3-20v  + 16v2)n2  + (l-8v 

m 

+ 8v2)y4  - 2(l-v) (3-12v  + 8v2)y2]  I I , > 
m m n-1  n+1 

+ ^ T’  { 2 [n3  - 2(5-10v  + 4v2)n2  + y2n  + 2 (1-v) (5-8v)n  - 4(1 

DR  mn  L m 

n 

-v)  (1-2v)y2]y2  I2  , + 2[-n3  - 2 (5-10v  + 4v2)n2  - y2n 
m m n-1  m 

- 2(l-v)(5-8v)n  - 4(1-v)(1-2v)y2]Y2  I2..  + 8[4(l-v)n4  + (3-10v 

m tn  n+1 

+ 4v2)Y2n2  - 4 (l-v)n2  + 2 (1-v) (l-2v)Y4]  I , I . , > 
m in  n-i  n+i 

+ 2^—1  ir-  {[-2n4  + 2 (2-3v)n3  - 3y2n2  - (3-32  + 8v2)n2  + 2(2 
mn  2 m 

DR  n 

? U i 7 L 122 

-v)y  n - y ]Y  I , + t-2n  - 2(2-3v)n  - 3y  n + (3-12v 

m m m n-1  m 

+ 8v2)n2  - 2(2-v)y2n  - yA]Y2  l2..  + 2[2n6  + (4y2  - 2)n4  + 3yV 
m mm  n+1  m m 

- ( 5+4 v - 8v2)Y2n2  + Y&]  I , I } 

tn  m n-1  n+1 


I 

I 

I 

1 

L 

l 

1 

1 

I 

1 


H.  - R2n“5  (l-v)n  {[-n2  - 4 (1-v) (l-2v)n  + 5-l2v  + 8v2jR2Y' 

in  Uq  J on 

+ [n3  + (l-2v)n2  - (l-2v)(3-4v)n  + 1 - 8v  + 8v2]RY  + 2(l-v)(n 

on 

-1)R2X'  + 2 (1-v) [-4vn2  - (l-4v)n  + 1]RX  - (n2-l)RT’ 

on  on  on 

+ n(n2-l)  T } 
on 

H2n  = 2^D~V~  R2n~3n  U~n2  " 2 (l-2v  + 4v2)n  + 3-12v  + 8v2]rVo 
o 

+ [n3  - (l-6v)n2  - (7-30v  + 24v2)n  - 3+12v  - 8v2]RY  - 2 [ (2 

on 

- 9v  + 8v2)n  - (1-v) (3-8v) ]R2x'  + 2[-4v(l-v)n2  + (2-v)(l-4v)n 

on 

+ (lv) (3-8v)]RX  + [-n2  + 2(l-4v)n  + (3-8v)]  RT'  + [n2 

on  on 

- 2(l-4v)n  - (3-8v)]n  TQn  } 

H3n  " r2TI~3  " (1_4v)n2  + 2 (1-v)  (~3+8v  - 8v2)n  + 2(1- 

o 

v) (3-12v  + 8v2)]R2Y'  + [-n4  + (l-4v)n3  + (l-2v)(7-16v  + 8v2)n 

on 

- 2(l-v) (3-16v  + 16v2)n  - 2(l-v)(3-l2v  + 8v3)]RY  + [(l-8v 

on 

+ 8v2)n2  + 8v(l-v) (l-2v)n  + 2 (1-v) (-3+12v  - 8v2)]R2X' 

on 

+ [8v(l-v)n3  + (-l-4v  + 24v2  - 16v3)n2  - 8v(l-v)n  + 2(l-v) 

(3-12v  + 8v2)]RX  + [n3  - (5-16v  + 8v2)n2  - 2v(3-4v)n 

on 

+ 2(1-vK3-8v)]RT'  + [-n4  + 2(l-3v)n3  + 2v(3-4v)n2  + (-3 

on 

+ 12v  - 8v2)]T  } 

on 

where 

D = ~J~3 — “ (2Y2n  [ (1-v)  (3-4v)n2  - v(5-4v)y2J  (I3^  + I3+1> 

R nJ 

+ y2  [n4  + (2y2  - 7+14v  - 8v2)n2  + , - I*,)  + 2n  [-8(1 

mm  m n— i n+i 


(24) 


-v)n  + (l+7v  - 4v  )y  n + 8(l-v)n  + v(5-4v)y  ] I , I 

m m n-l  n+l 

(I  . + I - ) + [-4n6  + (4-llY2)n4  + (-10y2  - 3+46v  - 40v2)y2n2 
n-l  n+l  mm  m 

- 3y6]I  . I (I  . - I } 

m n-l  n+l  n-l  n+l 

D = 2R3n"3  (l-v)(l-2v)(3-4v)(n2-l)n  (25) 

o 

The  symbols  Ir  and  In+^  appearing  here  represent  modified  Bessel  functions 

of  the  first  kind  of  argument  y = A R. 

m m 

The  particular  solutions  constituting  parts  of  eqs.  (21)  and  (22)  have  been 
found  by  applying  the  method  of  variation  of  the  constants: 

V1’"-"*  * - K™(I>  WV>  + Km‘r>  WV> 

- Tz3»n(r>  WV'  + ‘WV’ 

* TZWr)  Vl <V>  + Kn-l<V> 

*p(r,n,«,)  - Z5mn(t)  In(Amr)  - Z6nnl(r)  K^r) 


V (r,n,0)  - . . [r  Z (r)  - r 
rp  ’ ’ 4 (n+l)  Ion 


-(n+l) 


Z20„<r» 


+ 4(b)  (r”'1  Z3o„<r>  - Z4<,n<r» 


-(n+l) 


- 4(SI)-  <-r  Zlon<r>  + r ' ' Z2on(r” 


+ irbr  Ir"’1  z3c„<r>  - r'(n'1>  z4„„<r>J 


- *r-  Un  - r_n  Z,n_(r)] 


$ (r,n,0) 


(26) 


I 


Wr) 


• l |0ln(a)  - <-1)m  ' ^Jol  * <-D*  G,  <p)] 


Z2ma<r) 


/ IG  (o)  - (-!)■  G.  (o)  - G.  <p>  + (-1)"  G.  (p) ) 


[WV>/'WV)>d(’ 


. 


Z3mn(r) 


Z4,nn<r> 


/ (Gln(p)  - (-l)m  G2n(p)  + G3n(p)  - (-l)m  C4n(p)] 

lKn-l(Xmp)/Mn-l(AmP)]dp 

r 

/ [C.  <p)  - <-«”  G.(p)  + G,„  (p)  - <-l)m  C.  (p) J 


Wr> 


Z2p»“> 


Z3on(t> 


Z40„<r> 


7 [Gm(p)  " G2n(p>  - G3n(p)  + G4n(p)Jp‘n  dp 


/ tGln(p)  " G2n(p)  “ G3n(p)  + G4n(p)lpn+2  dp 


I !Gln(p)  - G2n(p)  + G3n(p)  - G4n(p)]p'n+1  dp 


7 [Gln(p)  - G2n(p)  + G3n(p)  - G4n(p)]pn  dp 


Z5mn(r) 


7 [Gln(p)  - (-l)m  G2n(p)][pKn(Amp)/Mn(Xmp)]dp 


Wr) 


7 IGln(P>  - (-1)"'  G2n(p)l[pIn(Xmp)/Mn(Xmp)ldp 


Z5on(r)  " 


J lG2n(p)  - Gln(p)]p_n+2  dp 
o 


I 

1 

I 

I 


and  where 


Z6on(r)  “ / tS»  - 


n+2 


2n 


In 


dp 


(27) 


M (X  p)  i K’  (X  p)  I (X  p)  - K (X  p)  I’  (X  p) 
nin  nm  n ta  nm  n m 


(28) 


As  it  is  seen  there  are  4 unknc  — functions  to  be  determined:  G^n(p),  •••*  G^n(p). 

To  this  end  one  has  to  substitute  the  expressions  (21)  and  (22)  to  the  in- 
version formulas  for  the  Fourier  transformation  with  respect  to  the  variable  z. 
Once  this  is  done  the  conditions  at  the  flat  ends  eqs.  13.1,  13.2,  14.1  and  14.2 
are  to  be  applied  (notice  that  condition  13.3  is  satisfied  automatically)  leading 
to  a system  of  Fredholm's  integral  equations  of  the  first  kind  with  respect  to 
aforementioned  functions.  Solution  of  this  system  will  provide  necessary  infor- 
mation for  the  evaluation  of  stresses. 
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APPENDIX  2:  Controlled  Impact  of  a Finite  Elastic  Rod  Upon  an  Isotropic, 
Elastic  Half-space. 

1.  Introduction 

The  present  paper  is  concerned  with  the  state  of  stress  in  a finite  elas- 
tic rod  striking  a stationary,  isotropic,  elastic  half-space.  The  motion  of 
the  far  end  of  the  rod  is  controlled  in  such  a way  that  the  depth  of  the  in- 
dentation is  within  the  limits  of  linear  elasticity.  In  order  to  avoid  the 
difficulty  resulting  from  the  fact  that  the  normal  displacement  of  an  elastic 
half-space  is  singular  at  the  point  of  application  of  the  concentrated  force 
it  is  assumed  that  the  reaction  of  the  rod  is  uniformly  distributed  over  the 
circular  area  corresponding  to  the  rod's  cross-section. 


2.  The  Solution  for  the  rod 

Let  u;  denote  axial  displacement  of  the  rod,  x - spatial  coordinate 
measured  along  its  axis  from  the  far  end,  and  E,  p - Young  modulus  and  mass 
per  unit  volume  of  the  rod's  material  respectively.  The  equation  of  motion 
of  the  rod  takes  the  following  form: 

O.xSL,  t>0  CD 


Let  the  function  f(t)  describing  controlled  motion  of  the  far  end  of  the  rod 
be  defined  in  such  a way  to  assure  that  the  displacement  be  zero  at  t = 0 

ui(0,t)  = f(t),  f CO)  = 0,  t > 0 (2) 

At  the  striking  end  of  the  rod  the  following  continuity  conditions  must  be 
satisfied: 


Ui(L,t)  = w(0,0,t)  = w(t),  0 * t 5 ti 
3ui  CL, t) 


(3) 


3x 


(A  +2y ) aw(0A-^  + A , 0 5 t 5 ti 

S 5 dZ  5 dr 


(4) 


I 
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where 

w(0,0,t)  = lim  w(r,z,t) 

u(0,0,t)  = lim  u(r,z,t)  for  r,  z ->  0 

are  the  vertical  and  radial  components  of  the  displacement  of  the  half-space 
at  the  impact  point  referred  to  polar  coordinate  system  r,  z.  The  axes  z 
and  x are  col  linear;  ti  denotes  the  separation  time  and  the  subscript  s des- 


cribes the  coefficients  related  to  half- space.  The  initial  conditions  for  the 
rod  are: 

/Hi,  IY  II  I 

x ^ L (5) 


Uj  (x, 0)  = 0,  — 1 = V = const.  0 < - < 


while  for  the  half- space  they  are 


u(r,z,0)  - w(r,z,0)  - » *(r,z,0>  . 0 


(6) 


The  elimination  of  the  variable  x from  eq.  1 is  achieved  by  applying  finite 

m = 1,  2,  ...  (7) 


Fourier  sine  transformation 

tl 


S{ui>=umW  = 

where  a = imr/L. 
m 


Ui (x,t)sin  am  x d x 


0 


The  transformed  eq.  (1)  takes  the  following  form 
d2um(t) 


nr 


* | ^ u ft)  = f a [f(t)-(-l)mw(t)] 
pmm  pm 


(8) 


0 < t < ti 

where  the  boundary  conditions  (2)  and  (3)  have  been  taken  into  account. 
Solution  of  the  eq.  (8)  with  consideration  of  the  transformed  conditions  (5) 
becomes : * 


um(t)  = V /p7E  = SE7p 


m 


sin  e (t-x)dx 
m 


[(f/T)  - (-l)mw(T)] 
0 < t < t 


0 


(9) 


where 


m 


am  /E7p. 


I 
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Applying  the  inversion  formula  one  obtains  from  (9) 

°°  1 - T 1 

Ui(x,t)  = (2VL/tt2)  /a/E  l — sin  e t sin  a x 

m=l  m mm 

rt  00 

+ (2/L)  /£/ 0 (f(x)  £ sin  e (t-x)  sin  a x 

j -1  hi  in 

0 m 1 

00 

- w(t)  l (-l)m  sin  em(t-x)  sin  am  x)  dx  0 < t 5 tj  (10) 

n*l 

Eq.  (10)  can  be  finally  put  in  the  form 

ui(x,t)  = (2VL/Tr2)A(x,t)  /p7E  + Ui(x,t)  - U2(x,t)  t < tj  (11) 

where  for  t 5 L/2  (t  = t /E7p) 

00  tt2x/2L  , x 5 t 

A(x,t)  = l — W-  sin  e t sin  a x = { TT2t/2L  , t < x * L - t (12) 

m=l  m mm  tt2  (L-x)/2L,  L-t  < x < L 

and  similar  results  can  be  obtained  for  t > L/2. 

Also  using  the  result 

<30 

Z sin  m X sin  mrrx/L  = j 5(x-^  x+LM)  for  M = 0,2,4,... 

and  similarly  for  M=l,3,5,...,  where  6 is  Dirac's  delta,  and  also  applying 
the  sifting  property  of  6 it  was  found  that: 

Ui(x,t)  = f (t  - x /p/E)  - f [t-(2L-x)  /p/E  ] 

+ f [t- (2L+x)  /p7E  ] - ...+(-l)M  f (t-  [(-1)M  x 

+ ML  + L]  /p7E  } t 5 ti  (13) 

U2(x,t)  - -w  [t-(L-x)  /p7E]  + w [t-  (L+x)  /p7E  ] 

- w [t-  (3L-x)  /p 7E  ] + . . . + (-l)M+1w  { t-  [ 

(-l)M+1x  ♦ (M+  )L  ]^7e  } t < ti  (14) 

The  displacement  (11)  depends  on  the  still  unknown  function  w(t). 


3.  Solution  for  elastic  half- space 

The  satisfaction  of  the  boundary  conditions  (3)  and  (4)  demands  that  the 
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vertical  displacement  w(t)  = w(0,0,t)  of  the  half-space  subjected  to  time 
varying  pressure  from  the  rod  be  determined.  Since  this  pressure  is  con- 
centrated at  a point  it  becomes  therefore  necessary  to  apply  SAITO  and 
CHONAN's  [1]  approximated  procedure.  Hence  it  is  assumed  that  the  rod's 
pressure  is  distributed  uniformly  over  the  disc  equal  to  the  cross-section 
of  the  rod.  The  contact  between  the  rod  and  the  half- space  is  maintained  -- 
as  originally  assumed  --  only  at  the  center  of  the  disc. 

As  an  auxiliary  problem  one  has  to  determine  vertical  deflection  wi  (r,z ,t) 
of  the  half- space  loaded  uniformly  over  a circular  disc  by  a suddenly  applied 
constant  unit  load.  The  solution  of  this  problem  has  been  obtained  by  EASON 
[2]  and  it  takes  the  following  form  for  r = z = 0: 


2 Try  .,  fn  n _ TrB2  , 02tt  / 1-a2  / a2-a2xz 

w,  (0,0, t)  - -grq-  ♦ — g s 


a 


- £ t y ✓ y2-l)(a2-y2x2)dy 

a { (i-^e2y2)2+  /Cy^Wj^TT 


£ | y(l-y2)  /(gV-l)  (a^-y2r2)dy 
a V0(l-^2y2)4  +(l-y2) (02y2-l) 


for  0 $ x 5 a 


(15) 


and  similar  expressions  valid  for  a 5 x 5 0a  etc. 

The  following  abbreviations  have  been  used  here 

02  3 (As  + 2ps)/us,  G = 202 (1-  ^02a2)  -/ 

- B2  , x = cit,  Cj  = /(Xs+2ys)/ps  (16) 

where  a is  the  coefficient  of  the  root  s = i aC  of  the  following  equation 

(€2  + ^02s2)2  - C2  /(r+rjU^s1)'-  0 (17) 

In  order  to  obtain  the  result  valid  for  a load  h(t)  varying  arbitrarily  in 
time  one  has  to  use  the  following  relationship 


(18) 


x(t)  = j\(0  x;  (t-c)d-  - jVt-ox;  (o  ^ 
o o 

relating  the  response  X(t)  to  arbitrary  excitation  h(t)  to  the  response 
Xj (t)  to  HEAVISIDE  type  of  excitation. 

The  resultant  force  at  the  striking  end  of  the  rod  is  obtained  from 
q(t)  = - E ^a2  t 5 t!  (19) 

where  eq.  11  should  be  substituted  for  ui(L,t). 

4.  The  conditional  equation 

Using  eq.  (18)  the  following  VOLTERRA's  integro-differential  equation 

with  respect  to  the  unknown  function  w(t)  is  obtained 

w(t)  = | q(t-OwJ  (0,0,0  dC  (20) 

0 

where  q(t)  depends  on  w(t)  through  eqs.  (19),  (11)  and  (14).  Integrating 

eq.  (20)  by  parts  one  gets  the  formula 

w(t)  = | wj  (0,0, t-C)  dC  (21) 

0 

The  LAPLACE  transform  of  (21)  takes  the  form 

w * p wjq  (22) 

where  p is  the  parameter  of  the  transformation  and  w is  the  transform  of 
w(t)  etc. 

After  numerous  rearrangements  the  following  explicit  equation  for  w 
results  from  eq.  (22): 

w (1  + /Ep  ir2a2p2Wi  [l+2e  2pb+2e  pb+...]  } 

* - pwiira2SEp  [*l+2e  pb-2e  2pb+...  ] 

♦ 2pT  [e’pb+e*3pb+...  ]}  (23) 


where  b * L /p/E. 


In  order  to  obtain  w(t)  what  is  necessary  for  the  analysis,  through 
eqs.  (11)  and  (14),  of  the  history  of  displacement  and  of  stress  in  the  rod 
one  must  invert  w resulting  from  eq.  (23).  This  can  be  achieved  by  using 
continued  fraction  representation  as  done  by  AKIN  and  COUNTS  in  [3]. 
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